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Abstract 



fi 

. The slave-boson mean-field method is used to study the two-dimensional 

', Hubbard model. A magnetic phase diagram allowing for paramagnetism, 

^ ■ weak- and strong ferromagnetism and antiferromagnetism is constructed and 

compared to the corresponding phase diagram using the Hartree-Fock ap- 
proximation (HFA). Magnetically ordered regions are reduced by a factor of 
^ ' awi ibout 3 along both the t/U and density axes compared to the HFA. Using 

■ the spin-rotation invariant formulation of the slave-boson method the helicity 

modulus is computed and for half-filling is found to practically coincide with 
that found using variational Monte Carlo calculations using the Gutzwiller 
wave function. Off half- filling the results can be used to compare with Quan- 
tum Monte Carlo calculations of the effective hopping parameter. Contrary 
to the case of half-filling, the slave-boson approach is seen to greatly improve 
the results of the HFA when off half-filling. 
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1 Introduction 



In the study of correlated electrons, for which both charge and spin degrees of free- 
dom are relevant, the Hubbard model is an intriguing simplification of reality that 
still contains a great deal of the essential physics Although superconductivity 
has not been demonstrated in this model, it is able to explain or reproduce a large 
number of experimental results on the copper oxides which superconduct at high 
temperatures 0. Despite such encouraging results, understanding of the model in 
dimensions two and higher is still rather limited. Even when the most elementary 
mean-field approximation, the Hartree-Fock approximation, is invoked, the phase 
diagram cannot be determined in full 0, since inhomogeneous phases, like spiral 
phases or domain walls, are able to supersede simple ferro- or antiferromagnetic 
phases. The possibility exists that more complicated phases not considered sofar 
are also important. Rigorous techniques like Quantum Monte Carlo and exact di- 
agonalization are limited to temperatures which may be too high and lattices which 
may be too small, respectively. Therefore it is of interest to employ approximations 
which go beyond the Hartree-Fock approximation (HFA). 

Some years ago Kotliar and Ruckenstein introduced, for the Hubbard model, 
the technique of using slave bosons to keep closer track of the site-occupancy than 
is done in the HFA. If a further approximation is made, the so-called "slave-boson 
mean-field" (SBMF) approximation, this approach was shown to be equivalent to 
the approximation scheme of Gutzwiller for the Hubbard model |Q. In a parallel 
development it was shown by Vollhardt and co-workers that the Gutzwiller ap- 
proximation scheme becomes exact in the limit of an infinite number of spatial 
dimensions, whereas the HFA does not become exact in this limit |^. Therefore, 
just like in classical statistical physics mean-field theory is a good starting point 
to study a specific model (because fluctuations become increasingly less relevant 
when increasing the dimension), for quantum models the SBMF approximation is 
a good starting point. In any case, it is expected to be an improvement over the 
HFA. In this connection, it is of interest to note that Oles and Zaanen compared 
the Gutzwiller approximation (GA) with the HFA in a two-band model of copper- 
oxide planes and showed that the GA is a good approximation for correlations 
at small length scales More recently the slave-boson approach of Kotliar and 
Ruckenstein has been refined in order to make it spin-rotation invariant [§, 0. 

In this paper, we make a detailed comparison of the SBMF and Hartree-Fock 
approximations to the (one-band) Hubbard model on a two-dimensional square 
lattice. Where possible we also compare to Quantum Monte Carlo calculations. By 
using the analytical result for the density of states of freely hopping electrons on the 
square lattice all calculations can conveniently be performed for a lattice of infinite 
size using one- dimensional integrals (over energy) only. Although we will mostly 
present results for the ground state (T = 0), the formulation we present is to a 
large extent valid for finite temperatures as well. A further advantage of the SBMF 
approach is that, in principle, it is valid for the whole range of Hubbard repulsion 
strengths and electron densities. First, we construct a ground-state phase-diagram 
allowing only for the simple magnetic phases: paramagnetic, antiferromagnetic, 
weakly and strongly (i.e. partially and fully polarized) ferromagnetic. We determine 
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all first-order and continuous transitions between these pliases.[] Although this 
approach may not give much information on the exact phase diagram (for instance, 
it is known that phases with spiraling magnetization supersede the antiferromagnet 
when going off half- filling) , it gives a clear impression of the improvement of the 
SBMF approximation over the HFA. Next, we derive an expression for the helicity 
modulus or spin stiffness in SBMF approximation using the spin-rotation invariant 
formulation and compare (for half-filling) to previous calculations of the helicity 
modulus in the HFA, as well as using the Gutzwiller wave function in a variational 
Monte Carlo calculation |TT|. It turns out that in SBMF approximation, apart 
from a negligible contribution, the helicity modulus is completely determined by 
the average kinetic energy (as is the case in the HFA) and therefore equivalent 
to the effective hopping parameter. When comparing the SBMF results for the 
effective hopping parameter to HFA and QMC calculations the improvement with 
respect to the HFA is only a few percent at half- filling and practically coincides 
with the results of variational Monte Carlo calculations using the Gutzwiller wave 
function, but the improvement is substantial off half-filling. 

The paper is organised as follows: in Section 2, we briefly introduce the slave- 
boson mean-field method for the Hubbard model and present the free energy for 
the antiferromagnetic, ferromagnetic and spiral phases as well as the corresponding 
consistency equations. In Section 3, we construct the ground-state phase diagram 
for the simple magnetic phases listed above and compare to the corresponding HFA 
phase diagram. An expression for the helicity modulus is derived and its connection 
to the effective hopping parameter is discussed in Section 4. The SBMF results are 
compared to both the HFA and QMC calculations where possible. The last section 
contains a discussion of the results and draws some conclusions. 



2 Slave Boson Mean Field method 

The Hamiltonian for the Hubbard model is given by: 

ij(T i i(7 

where creates an electron at site i with spin a, = cl^Ci„, tij is the one- 
electron transfer integral between sites j and i {tij equals t if i and j are nearest 
neighbours and otherwise), U the on-site repulsion [U > 0), and fi the chemical 
potential (/i = U/2 corresponds to a half-filled lattice, i.e. n = J^iai'i^'ia) = !)• 
The Hubbard model allows for four different occupancies of a single site: it can 
be empty, singly occupied by either a spin-up or spin-down electron, or doubly 
occupied (by electrons of opposite spin). This leads to the idea of introducing 
four different kinds of "slave" bosons, one for each of the possible occupancies. In 
order for this to become a bookkeeping device one introduces the constraints that 
at each site there is always exactly one boson present and that it is of the kind 
corresponding to the electron occupancy. If e, pi and d denote the annihilation 

^Some of the results on the SBMF phase diagram were also obtained in Ref . pOj . 
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operators for the four kinds of bosons, these constraints are: 



ejci + Ea p\aPia + = 1 for all i, (2) 
claCia = p\aPia + d\di for all 1 and a =t, j, (3) 

The interaction term in the Hamiltonian (|ID can now be replaced by one containing 
only the counting operator for (i-bosons. In order for the boson presence to keep 
in correspondence with the electron occupancy for each site, the hopping term in 
the Hamiltonian needs to be adjusted by connecting to each electron annihilation 
operator Cjo- the following boson-transformation operator: 

= elpia + pl_^di. (4) 

In fact this choice for z is not unique and we follow Ref . |Q by making a choice 
which in the case of [/ = gives the correct result if a subsequent saddle-point 
approximation is made: 



la 



1 - dldi - ptpia) ' Zia (l - eld - pl^^PL-a) ' • (5) 



In the physical subspace, defined by of the enlarged boson-electron Hilbert 

space, the Hamiltonian, 

Ti-SB = - 51 UjC^ztzjaCja + UY^ d\di - yu X! ^^'^^ (6) 

ijcF i i(T 

has the same matrix elements as the original Hamiltonian in the original Hilbert 
space containing only electron states. Therefore, up to here only a reformulation 
of the original problem has been achieved. However, if now in the functional inte- 
gral formulation the saddle-point approximation of time- and position-independent 
Bose fields is made, a set of equations results which is similar to those found in the 
Hartree-Fock approximation, but more general. This approximation is called the 
Slave-Boson Mean-Field (SBMF) approximation. For comparison, the HFA can 
be obtained in such a functional integral formulation by first applying a suitable 
Hubbard-Stratonovich transformation and then making the saddle-point approxi- 
mation for the original problem For further details on the functional integral 
formulation for the Hubbard model and the subsequent saddle-point approximation 
we refer to previous papers on this subject 0, ||, ||, |14|, p!6| . 

In the next three subsections, we present the expressions for the free energy 
of the ferromagnetic, antiferro magnetic, and spiral phases. For the ferromagnet 
and antiferromagnet we also give the consistency equations for the appearing self- 
consistent fields, which must be obeyed for the free energy to become minimal. To 
some extent this is a repetition of previously published results [rU, 0, but they 



are given for the reader's convenience and to establish the notation. 



2.1 Ferromagnetic Phase 

For the ferromagnetic phase one assumes a non-zero homogeneous magnetization m, 
even if a magnetic field is absent. To be able to calculate the magnetic susceptibility 
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we also include a magnetic field h in the Hamiltonians Ti. and TisB by adding a term: 



n 



mag 



(7) 



Here and in the following we adopt the convention that if a does not appear as 
an index it attains the values +1 and —1 if the corresponding index is j and |, 
respectively. In order to treat the density of electrons n and magnetization m on 
an equal footing, we introduce a slightly different definition for the free energy per 
site, to be denoted by ip. For the ferromagnet our free energy is defined as: 



(3N 



lnTre"'^^'^s''+^'""'5^ + Hn + hm . 



It is given by (see also Refs. 

1 



Ud + fin + Am 



where 



with 



k,(T 



EM 



(8) 



(9) 



(10) 



t(k) = — 2t(cos kx + cos ky) , 

the band structure of freely hopping electrons on a square lattice. The density n 
and magnetization m are given by: n = + ni and m = rij — n| with rio- = (rijo-)- 
Since in principle the electrons have been integrated out in obtaining (P), n and 
m can be understood as a shorthand for the Bose fields p| and p| via the relation: 
rio- = + (which is the average of constraint (H)). The parameters Ji and A 
are an effective chemical potential and effective magnetic field, respectively, which 
incorporate the Lagrange multipliers A^^-* used to enforce the constraints (^): 



A 



h 



{\f + A 



A 



(12) 
(13) 



The Lagrange multiplier A^^^ associated with constraint (0) has disappeared because 
the average of this constraint must hold. In fact, both constraints are only satisfied 
on average in the saddle-point approximation. The band-renormalization factor go- 
appearing in (|10D is in this approximation of time- and position-independent Bose 
fields a function of n,m and d (as follows directly from (|)-(||)): 



q^{n,m,d) = {z]^Zi„) 



1 — n + d'^){n + am — 2d'^) + d^/n — am — 2d'^ 



{n + am) [1 — i(n + am)] 



(14) 

We note that if this is rewritten as a function of n^, and d one exactly recovers 
the expression for the band renormalization in the Gutzwiller approximation (see 
Ref . jlTl ; our d"^ is called d there0). This would not be true if another choice for Zjg- 
than d^) had been made. 

^In Ref.M similar renormalization factors are derived in an alternative manner. 
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The sum over k in is over the whole Brillouin zone of the square lattice. In 
the limit of an infinitely large lattice, which can be treated after the approxima- 
tions made, the resulting integral over the (two-dimensional) Brillouin zone can be 
rewritten as an (one-dimensional) integral over energy using the density of states 
(DOS) of freely hopping electrons. For the square lattice this DOS, Af{e), is known 
analytically: 



k 







\e\ < At 
\e\ > At 



(15) 



where K{x) is the complete elliptic integral of the first kind [|n 
DOS the free energy is: 



Employing the 



--E 



deATie) In 



1 + e 



+ UdP + nn + Am 



(16) 



with Efj{e) = q^e — aX — ft. Unless stated otherwise explicitly integrals over e run 
from — oo to oo (the relevant integration range is of course limited by (|15|)). 

For a given interaction strength U, ifp is now given as a function of the five 
variables n,m,d,fi and A, whereas /i and h are control parameters, regulating (al- 
though not directly in this slave-boson approach) the density and magnetization. 
The optimal values for the five variables must be found from the three minimiza- 
tion conditions: dip^/dd = 0, dipp/dfi = 0, dipp/dX = 0, as well as from the two 
equations arising from the Legendre transform between grand potential (function 
of fi and h) and free energy (function of n and m): dipp/dm = h and dipp/dn = /i. 
Applying these conditions to (|16|) results in: 



U 
n 
m 
X 



2d ? ^'^"^'^ 

cr 

(7 

1^ ^ ^ Qcrn^a j 



(17) 
(18) 
(19) 
(20) 
(21) 



where go-o denotes the first partial derivative of with respect to a (= n, m, d) 
and we introduced the abbreviations: 



n„ = I deAf{e)f[E^is)] , 



(22) 
(23) 



with the Fermi-Dirac distribution, f{E) 



1 + e' 



The five equations (p^)- 



(pTD need to be solved self-consistently for given ?7, jj, and h and the results can be 
inserted in (0) to obtain the corresponding free energy. In practice, we will often be 
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interested in calculations for a fixed density n, in which case the last equation (pT]) 
does not appear (/i disappears from the problem, however fi remains). Expressions 
for partial derivatives of q„ are given in Appendix A. 

The above allows to compute the free energy of three different phases: the 
paramagnetic phase (PM), for which m = if /i = 0, the strong ferromagnet (SF), 
for which m = n and both are non-zero even ii h = 0, and the weak ferromagnet 
(WF), for which m < n and both are non-zero even if /i = 0. In Section 3, we 
will obtain the lines in the {t/U, ?T,)-diagram of first-order phase transitions between 
these phases. In order to obtain the line of the continuous phase transition between 
paramagnet and ferromagnet (the ferromagnet can be either strong or weak), one 
needs to find where the susceptibility x of the paramagnet diverges. An expression 
for X is derived in Appendix B. 



2.2 Anti-ferromagnetic Phase 

For the antiferromagnetic phase we divide the square lattice in two sublattices, such 
that points on one sublattice have only points of the other sublattice as nearest 
neighbours. Furthermore, we assume a non-zero staggered magnetization rris, i.e. 
the magnetization is rris on one sublattice and — mg on the other. To be able to 
calculate the staggered susceptibility we add a staggered-magnetic-field term to the 
Hamiltonians 7i and ?^sb: 

where hi s equals on one sublattice and —h^ on the other. The saddle-point 
approximation of time- and position-independent Bose-fields on each of the two 
sublattices separately (introducing staggered Lagrange multiphers as well) results 
in a 2 X 2 problem (for each a separately, only one-electron states with k and 
k-|- Q couple, Q = (tt, tt)) for the quasi-particle spectrum, which is easily diagonal- 
ized. Analogously to the ferromagnet, we define the "free energy" per site for the 
antiferromagnet as: 

^ In Tre"^^^ss+^'"^S'=) + ^ln + h^m^ . (25) 



(3N 



^ ' In [l + e-^^e^)] + Ud^ + fm + \,m, , (26) 



V^AF = - 

It is given by: 
where 

Eiy) = ±4^iim + k - . (27) 

and we have again introduced an effective chemical potential fi and an effective 
magnetic field Ag. The prime indicates that the sum over k is over the magnetic 
Brillouin zone only [k^ i: ky E [— tt, tt]). The band-renormalization factor Qs is now 
a cr-independent quantity because of the staggering: 

gs(n, ms, d) = {zj^Zi^) af = ZAaZsa = zi^z^ , (28) 
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where 

:i-n + d^) in + am,-2d^) + dVn-am.-2d'^ ^^^^ 

n + arris) (l - ^^±f^) 

As for the ferromagnet, the free energy can be expressed as an integral over the 
DOS of freely hopping electrons, Af{e): 

(PAF = ~J deAfie) In [l + e'^^^")] + Ud^ + fin + ~X,m, . (30) 

For convenience we restrict ourselves to densities n < 1; in that case only the 
negative square root in (^) is relevant and we have in 



Eie) = -^q!e^ + X^^~jl. (31) 

The consistency equations for the antiferromagnet are obtained in the same way as 
for the ferromagnet and read: 

f/ ^ g., (32) 
n = j deM{e)f[E{e)], (33) 

rris = As / de ) ^ , (34) 



where we have defined: 



'q^e^ + A2 

As = hs + qsm,£ , (35) 
fi = IJ^ + qsn£ , (36) 



,,^($£fm. (37) 



Expressions for partial derivatives qsa of q, are given in Appendix C. (Paf is a 
function of the five variables n,ms,d,fl and As, whereas fi and are control pa- 
rameters. Calculating (p^F for fixed U/t and n, means that the last consistency 



equation (36) becomes irrelevant again and the four remaining variables must be 
found self-consistently from the four remaining consistency equations. To compute 
the free energy the staggered field hs is taken to be zero. This free energy may be 
compared to the free energies of the three phases discussed in the previous sub- 
section and lines of first-order transitions in the (t/[/, n) -plane may be found (see 
Section 3). We note that the equations for the paramagnet can also be found from 
the above equations for the antiferromagnet by putting rris = if /is = (then 
also As = 0). In order to find the transition line for the continuous phase tran- 
sition between antiferromagnet and paramagnet an expression for the staggered 
susceptibility Xs is derived in Appendix D. 
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2.3 Spiral Phase 



To obtain the spiral phase the magnetization vector is assumed to vary in space as: 

nij = m (cos(q • Rj), sin(q • Rj), 0) . (38) 



In Ref . [p^m , the spin- rotation invariant formulation is used to compute the free 
energy ipsp for this phase. For details of this calculation we refer to Refs.p], |16 
here we suffice by quoting the result (in our notation): 



where 



V^sp 



(3N 



l + e 



-/3Eq,.(k) 



+ Uii^ + iln- A(2) 



m . 



(39) 



f(k)+t(k+q) 
2 



/i + 



f(k)-t(k+q) 
2 



+ 



_[t(k)+t(k+q)]+A(2) 



1 

21 2 



(40) 



with z/ = ±1. The parameters z± are functions of n,m and d; in terms of the z^ 
(formula ([29| ) with rris replaced by m) they are given by: 



z± = I (^T ± ^i) • 



(41) 



The parameter fi = fi — X'^^ is again an effective chemical potential, whereas Aq^"* 
and A^^-* are the Lagrange multipliers arising from the constraint (|^) when made 
spin-rotation invariant .0 We did not include any explicit magnetic field in the 
Hamiltonian in studying the spiral phases. One may verify that for q = (0, 0) and 
q = (tt, tt) ( ^OD reduces to the expressions ( |10|) and (|2^) for the ferromagnet and 
antiferromagnet, respectively. In terms of the z± the band-renormalization factors 



are given by Qa 



+ az^) and Qs 



zl — zz 



For the spiral phase we do not give the consistency equations like we did for the 
ferromagnet and antiferromagnet, since we are only aiming at a simple phase dia- 
gram which does not include the spiral phases. Moreover, the consistency equations 
and free energy cannot be expressed as one-dimensional integrals over a density of 
states because of the spiraling vector q involved. Therefore, the consistency equa- 
tions need to be solved numerically on a finite lattice. The regions in the phase 
diagram where spiral phases dominate the simple magnetic phases were calculated 
in Refs.|j6[ and [|19|. In this paper, we will only use the expressions above to derive 
a formula within SBMF approximation for the helicity modulus (or: spin stiffness) 
and effective hopping parameter in Section 4. 



3 Magnetic Phase Diagram 

We compute, in SBMF approximation, the complete (i.e. all first order and contin- 
uous phase transitions are included) ground-state magnetic phase diagram for the 

■^In the spin-rotation invariant formulation, constraint ^ gives rise to a scalar Lagrange mul- 
tiplier Aq^'' as well as a vector Lagrange multiplier A'^^. For a spiral phase the latter results in 
another scalar A^^^^ = |A'^^-'| because it must show the same spatial variation as m in (^). 
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Hubbard model on a square lattice allowing for the four simple magnetic phases, 
paramagnet (PM), weak ferromagnet (WF), strong ferromagnet (SF) and antifer- 
romagnet (AF). In their original paper, Kotliar and Ruckenstein ^ only calculated 
the lines where the PM becomes unstable towards ferromagnetic or antiferromag- 
netic ordering (continuous transitions), whereas Evans [^] also included some first- 
order transitions, but not all, so that an incomplete picture emerged. 

First, the regions in the (4t/f/, n)-plane where ferromagnet ism (either WF or SF) 
and antiferromagnetism can occur are determined, by calculating the lines where 
the homogeneous and staggered susceptibilities, x and Xs^ of the PM diverge. In 
Appendices B and D expressions for x ^"^^ Xs are derived. The condition that the 
denominator in these expressions vanishes (generalized Stoner criterion) provides 
an additional equation to be solved in conjunction with the consistency equations 
for the PM (see Section 2.1). In this way, for fixed n, the additional equation fixes 
the (critical) U /t value for which the susceptibility diverges. The resulting lines are 
displayed in Figure 1(a) and agree with previously published results |10 . 



Now, using the formulae in Sections 2.1 and 2.2 all first order phase transition 
lines in the (4t/?7, n) -diagram are computed. Since no susceptibilities are required, 
h and /is are taken to be zero. In principle, for each of the four phases for fixed 
values of U/t and n the energy is found by solving the consistency equations si- 
mult aneously.^ For the SF and PM this problem simplifies somewhat: for the SF 
the set of equations (|l^)-(^) is reduced by one (since m = n) and for the PM we 
have m = 0. For each pair of phases one then finds a line in the (4t/?7, n)-plane 
where the two energies are equal. The results of such calculations are displayed in 
Fig. 1(a). In principle there are six such lines, but the first order PM/AF transi- 
tion line coincides with the continuous PM/AF transition. Note however that the 
continuous PM/F transition and the first order PM/WF transition differ. 

Taking into account all four phases, the phase diagram of the Hubbard model on 
a square lattice in SBMF approximation of Figure 1(b) emerges, in which all inter- 
rupted lines denote first-order transitions and the full line a continuous (PM/AF) 
transition. We now discuss the phase diagram in comparison with the same phase 
diagram as obtained in the Hartree-Fock approximation (HFA) and in comparison 
with previously published SBMF results. 

The corresponding, i.e., allowing for the same four phases, phase diagram to 
Fig. 1(b) in the HFA is shown in Figure 1(c). A similar diagram was given previously 
in three dimensions by Penn [21] and in two dimensions by Hirsch [2^, but in the 
latter the non- monotonous behaviour of the F/AF transition line was missed and 
the region of WF was not determined. Long has given the PM/F/AF phase diagram 
in HFA using a constant density of states; in that case no extremum in the F/AF 
transition occurs [^. Another surprising feature of Fig. 1(c), besides the maximum 
in the F/AF line, is the fact that the WF/SF transition line is found to oscillate 
slightly around the line of the continuous P/F line. The difference in these curves 
is very small, but we have ascertained that it is not due to numerical inaccuracies. 
Thermodynamically such behaviour is allowed; it only means that along the P/F 

""^Using the analytically known density of states (jlj) this may be conveniently done with the 
program MATHEMATICA [Q. Some care is required in integrating through the logarithmic 
singularity of A/'(e) in e = 0. 
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boundary the transition is sometimes continuous and sometimes first order. In 
comparing Fig. 1(b) and Fig. 1(c) it should be noted that they are topologically the 
same, but that because of the difference in scale on both the density and At/U axes, 
the SBMF has reduced the magnetically ordered regions considerably with respect 
to the HFA. Furthermore, the region where WF dominates has grown at the expense 
of the SF phase. In Table I, we list the location of the "tripod" points (i.e., points 
where three phases meet; the AF/SF/WF and PM/SF/WF points are triple points, 
i.e., points where three first-order transitions meet), as well as the renormalization 
factor obtained in going from the HFA to the SBMF approximation. Globally 
speaking this factor is about 3 for the hole density 5 and also about 3 for t/U (if 
the ferromagnetic region is considered as a whole). Since the HFA overestimates 
the importance of magnetic ordering p3|], the SBMF approximation is clearly an 
improvement. The critical hole density above which antiferromagnetism cannot 
occur is determined by the continuous PM/AF transition and is given by 5^^ = 0.21 
in SBMF. More interestingly, the critical hole density above which ferromagnetism 
cannot occur is in SBMF determined by the first order PM/SF transition and 



given by exactly 6^ = 1/3, as a simple argument can show (see e.g. Ref. ||19||). The 



latter value agrees very well with the result 5f = 0.29 obtained from calculations 



using a variational wave function |24|. Remarkably, also high-temperature series 



expansions for the Hubbard model find for f//t — > oo a value of about 0.33, below 



which ferromagnetic nearest-neighbour correlations occur [25|, and a value "near 



3/11" (~ 0.27), below which a strong separation of energy scales for spin and 
translational degrees of freedom is observed |^. Although these features in the 
high-temperature series expansions appear to be temperature independent over a 
wide temperature range, extrapolation to T = is cumbersome in such expansions 



2^ . We also note that the x~ = line is nowhere in the diagram a phase boundary; 
therefore the continuous PM/F transition, in this approximation and contrary to 
the HFA result, is preempted by first order PM/WF or PM/SF transitions. 

A calculation of the phase diagram similar to ours was previously performed 
by Evans [IC]. However, the more cumbersome PM/WF and AF/WF first order 
transitions were not computed and for the WF/SF transition only the limit m = n 
and (i = in the WF was taken. The latter determination turns out only to give an 
upper bound (in t/U , for fixed n) for the WF/SF first-order transition computed by 
comparing energies, as it should. As a result the final phase diagram of Ref. is 
obtained by removing from Fig. 1(a) the PM/WF and AF/WF lines and replacing 
the WF/SF line by one extending from (0.15,0.0) to (0.0,0.38) in the (4t/f/,(5)- 
plane (5 = 1 — n) . Evans then appears to call WF only the tiny, triangle-shaped, 
region at the center of our WF region (although this is not very clear from Fig.l 
in Ref. |]10[)- This assignment, however, is thermodynamically not justified since 
two of the boundaries then correspond to PM/SF and PM/AF transitions. Also 
the third boundary in that case (PM/F) is not a boundary for the WF region as 
computed by us. 

To conclude the discussion of the phase diagram in Fig.l, we stress that the 
actual ground-state phase diagram of the Hubbard model, even when constructed 
within either the HF or SBMF approximation, will also have to include inhomo- 
geneous phases like domain walls and spiral phases. For instance, it was shown 
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within the SBMF approximation that spiral phases supersede the antiferromagnet 
immediately when going off half-filling and also the ferromagnetic phases shrink 



somewhat in favour of certain spiral phases [0, 19|. Our detailed determination 
of the simple phase diagram only serves to study the consequences of the SBMF 
approximation when compared with the HFA. 



4 Helicity Modulus and Effective Hopping 

A crucial quantity in the study of quantum-ordered states is the helicity modulus, 
which is the stiffness associated with a twist of the order parameter, or, equivalently, 
with phase fluctuations of a complex order parameter. For the attractive Hubbard 
model, which exhibits superconducting or superfluid order, the helicity modulus 
corresponds to the superfluid density, whereas for the repulsive Hubbard model it 
is the spin stiffness of the AF ordered phase at half-fllling. In previous papers, the 
helicity modulus, denoted by ps, was calculated for the 2D Hubbard model both in 



the HFA as well as by variational Monte Carlo methods |TT], [T^. Also a comparison 
with exact diagonalization and Quantum Monte Carlo results was made, showing 
that the HFA renders quantitatively reasonable results for ps |jl3|. For the repulsive 



Hubbard model this could only be shown at half-fllling. Since we have already seen 
that the SBMF approximation is an improvement over the HFA, it is of interest 
to see what it will give for p^. In this section, we flrst derive an expression for p^ 
within the SBMF approximation and use it to compute ps- Then the connection 
between ps and the effective hopping parameter is discussed, leading to calculations 
of the effective hopping both at half-fllling and off half-fllling. The SBMF results 
are compared to results from the HFA and Quantum Monte Carlo results. 

To obtain ps, we can make use of the results for the spiral phase in Section 2.3. 
In the AF phase the order parameter is the staggered magnetization. This can be 
viewed as a spiral phase with spiral vector q = Q = (vr, vr). A small twist in the AF 
order parameter then corresponds to a spiral vector which deviates slightly from 
Q: 

q=(7r,7r)-q. (42) 
The helicity modulus ps is given by: 



Ps 



lim 



ig2 



(43) 



with the free energy per site (p given by ( pOf ) and where q is the modulus of q. To 
facilitate the computation it is advantageous to perform a further manipulation. 
Since the spectrum is periodic in reciprocal space and we are going to integrate 
over the full Brillouin zone, it is allowed to shift the spectrum over iq, so that ( ^OD 
becomes (if in turn we rename iq, for notational convenience, to q): 



f(k+q)-f(k-q) 
2 

t(k+q)+f(k-q) 
2 



P + 



+ 



z+z_ [t(k + q) - t(k - q)] + A(2) ,(44) 



1 

2^ 2 
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with V = ±1. After this manipulation, the occurring sum t(k + q) +t(k — q) is even 
in the small parameter q (=|q|) and the difference is odd. The same is not true for 
the small parameter q in the occurring sum and difference t(k) ± t(k + Q — q) in 



(^Ol). In this notation, is given by: 



P^ = }^o ^2 ' (45) 

If we now restrict to ground-state properties (T = 0) and densities less than half- 
filling (n < 1), we only need to expand the u = — 1-branch for small q (and integrate 
over the Brillouin zone (BZ)) to obtain (f{q) — f{0) to order q"^. We find (taking 
q = (g, 0) for convenience) that the term proportional to q vanishes after integrating 
over BZ, and that ps is given by: 



2 



t(k)cos(A;^) AtzlzH"^ {lijsm'^ih 



+ 



2E{k) E^k) 



(46) 



with 



Eik) = ^g2^2(k) + X2 , (47) 
and the band renormalization q^ is: 

qs = zl-z^_. (48) 

We have omitted the superscript (2) on A. We remark that if the above shift in 
the BZ is not performed, a much longer expression for p^ results; the expression is 
equivalent to (^), but this is not trivial. We further note that the T = HFA result 
of Ref.[0 is recovered by omitting the second term in the BZ sum and putting qs 
equal to 1. 

In order to compute ps for fixed density n, according to Section 2.2, the following 
set of equations needs to be solved self-consistently (for T = and = 0; cf. 

my. 

m, = 2A /'d£— ^^^feL , (49) 
A ^^2 + A2^ 



- 2 /-^ Af(e)e^ 

A = -Qsm, / de —j^ , (50) 

^ = -f ;T7^ ' 51 



(^2 + A 

where A = A/gs and /2 is determined by the fixed n: 

n = 2 deAf{e) , (52) 

J u 



In terms of the parameters of Section 2.2, we have p = \l p? — Xl/q^i (As is called A 
here). The band renormalization q^ is given by (p8D-(p9D. In terms of an integral 
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over the DOS, the energy of the AF state (per site) and the spin stiffness ps are 
given by: 



CAF = -2gs r d£ -^(^) \/^' + + Ud^ + \q,m, (53) 

J fl 



2 



P^ = J - -TT72 - - ^ ,3/2 (54) 



where A/'v(£:) is the weighted density of states: 



-A/'v(^) = ^E[Vt(k)]^5(e-t(k)), (55) 



which for a square lattice can be calculated analytically (see Ref.[l28| 



AC 




1 , 

At 



At 



K 



(56) 



for \e\ < At and zero otherwise. K{x) and i?(x) are the complete elliptic integrals of 
the first and second kind, respectively. We remark that using the weighted density 



of states the finite-temperature result for ps obtained in the HFA can also be 
written as an integral over energy. 

Since only at half-filling the AF phase is the ground-state, we first restrict 
ourselves to this case: n = 1 (p = 0). In Table II, self-consistent parameters are 
given for various U /t\ the corresponding energy caf and band renormalization 
are also given. These results agree with those published previously by Hasegawa 
[ p9[| . In particular, we note that q^ never deviates from 1 by more than 5 %, 
implying that slave bosons at half-filling renormalize the Hartree-Fock results only 
by a small amount. As concerns p^-, for the case of half- filling the integral containing 



M^{e) plays no role since for n = 1 we have Z- = Q (as is easily verified from (^9]) 



and (0)). In Table III and Figure 2, we compare the results for obtained in 
the SBMF, with those obtained previously using the HFA and using variational 
Monte Carlo calculations with an (antiferromagnetic) Gutzwiller wave function 
(GWVMC) The fact that slave bosons only renormalize the HF results a little 
bit is reflected in the fact that the SBMF and HF results never differ by more 
than 7 %. Note however that the SBMF result is always larger (except for very 
small f//t, U/t < 2.5), whereas q^., which enters as a factor in (|5^) , is smaller than 
1. The direct effect of q^ is more than compensated by a renormalized (smaller) 
value of the antiferromagnetic gap A = Ags- A further observation from Table HI 
and Fig. 2 is that the SBMF results almost coincide with the GWVMC results. 
The difference between the results is an indication of the difference between the 
Gutzwiller approximation (which is equivalent to the present SBMF approximation, 
see the introduction) and the Gutzwiller wave function. Although these are not 
identical [^], the difference for the spin stiffness is not big, as shown in Fig.2. 
Therefore the tedious variational Monte Carlo calculations for ps can be replaced 
by the above set of equations which are exact (within the SBMF) and easy to solve. 

On very general grounds it can be derived that the helicity modulus (spin stiff- 
ness for positive U and superfiuid weight for negative U) comprises of a "direct" part 
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proportional to the average kinetic energy (T) and a part related to the current- 
current correlation function A^^. The HFA effectively neglects the A2.j.-part, 



whereas GWVMC calculations find only a negligible correction to the kinetic part 
— i(T) 0]. Since in formula (|46D the first term is exactly — |(T) in the SBMF 
approximation and at half-filling the second term equals zero, we can conclude 



that also the SBMF approximation only gives the kinetic part of ps- In Ref. ||13|| 
it was estimated (by comparing to appropriate exact diagonalization calculations 
and Quantum Monte Carlo (QMC) calculations) that the HFA overestimates ps at 
half-filling by 68, 40, 38 and 36 % for U/t = 4, 8, 10 and 20, respectively. However, 
if one compares the kinetic energy found in the HF and SBMF approximations with 
QMC calculations for n = 1 (The latter are obtained from Ref.[^), one concludes 
that the approximations perform very satisfactorily. This is illustrated in Figure 3, 
where we plot the effective hopping integral teff, defined by normalizing the average 
kinetic energy for interaction constant U with that for U = 0: 

teS _ {(^ia^jcr + Cj^Cia)u (57) 

The denominator is easily evaluated as the energy of the PM phase in the HFA, since 
for U = the HFA is exact and there is no potential energy. The QMC data is taken 
at sufficiently low temperature (/3t = 16) for this comparison with T = results 
to be meaningful. We note that a similar comparison of SBMF and QMC data 



was made in Ref.[^; in that paper the SBMF approach was formulated as a 14- 
dimensional optimization problem. By comparing their Figure 2(b) with our Fig.3, 
the results from our simple formula (|5^) are found to be the same. A surprising 
feature of Fig.3 is perhaps not so much that the SBMF results approximate the 
QMC results so well, but that the HFA results do the same already. 

We now discuss the off-half-filling case; since off half-filling their exists a spi- 
raling vector qo for which the spiral phase has lower energy than the AF phase, 
we cannot call ps as given by (0) the stiffness of the ground state anymore. In- 
stead, the expression ( ^4|) has the interpretation of the stiffness of the AF phase 
with respect to a small deviation from spiraling vector (7r,7r). In order to com- 
pute the stiffness of the ground state one would have to perturb the spiral phase 
with vector qo, but this is beyond the scope of the present paper. Here we only 
investigate how well the effective hopping (or, equivalently, the kinetic energy) of 
the two-dimensional Hubbard model off half-filling is described by the first term 
in (^4]). In Figure 4(a), we compare tes/t obtained from Hartree-Fock and SBMF 



approximations (for T = 0) with low-temperature {(3t = 6) QMC data (The latter 
are obtained from Ref.[^). The results are displayed as a function of density n 
for the one value of U/t (=4) for which there are QMC data available. From the 
phase diagrams in Fig.l it is clear that both in the HFA and SBMF approximation 
an AF/PM transition occurs^ for some critical density ric [ric = 0.86 in SBMF and 
He = 0.76 in HFA for U/t = 4, the former is not showing in Fig.l(a)). Below ric 
(paramagnetic phase), tes/t equals 1 in the HFA and equals the band renormal- 
ization q {=qi = Qi) in the SBMF approximation. Clearly the SBMF approach 

^If the restricted set of four phases is considered as before. 
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is a significant improvement over the HFA; the agreement with the QMC data is 
less good in the density interval just off half-filling. This is most probably caused 
by the fact that for these densities the assumed (antiferromagnetic) phase in the 
SBMF approach is not the correct one. The same remark concerning Ref. |3T| as 



made above for n = 1 is appropriate here. Finally, in Figure 4(b), we also show the 
results for the effective hopping integral for a few other values of f//t, for which no 
QMC data are available. The corresponding HFA curves are not shown, but the be- 
haviour is similar as for U/t = 4: at half-filling tcs/t is somewhat below the SBMF 
result and rises to 1 in going off half- filling. The densities below which t^g / 1 equals 
1 can be read off from Fig. 1(c). We note that in Fig. 4(b) the non-differentiability 
at Tie (which is close to 0.8 for U/t = 8, 12, 16, as can be seen in Fig.l(a) and (b)) 
is less pronounced for [//t = 8, 12, 16 than it is for U/t = 4. 



5 Discussion and Conclusions 

Above we have given a detailed account of calculations within the slave-boson 
mean-field (SBMF) approximation for the repulsive Hubbard model on a square 
lattice. We have focused on the phase diagram, a particular response function, the 
helicity modulus ps, and the related effective hopping integral. All calculations can 
be expressed in terms of a set of integral equations with one-dimensional integrals 
over energy containing a density of states, which is known analytically for the square 
lattice. These equations are solved self-consistently. 

If for the {t/U,n) phase diagram we restrict to simple magnetic phases, the 
SBMF approach is found to reduce the magnetically ordered regions with respect 
to the Hartree-Fock approximation (HFA). Along the density axis the reduction is 
roughly a factor of 3, whereas along the t/[/-axis the reduction of the ferromag- 
netic region (weak- and strong ferromagnetism together) is also a factor of 3. In the 
SBMF approach the portion of weak ferromagnetism grows at the expense of the 
strong-ferromagnetism portion when compared to the HFA. The present SBMF 
phase diagram is more likely to be a good starting point for more sophisticated 
approaches to the phase diagram than is the HFA phase diagram, because in an 
infinite number of dimensions the SBMF approach becomes exact (see the intro- 
duction). On the other hand, it is seen that the SBMF approach is not qualitatively 
different from the HFA, but rather a renormalized form of HFA. 

The new quantity that we obtain within the SBMF approach is the helicity 
modulus Ps. At half-filling, the results for ps practically coincide with those obtained 
using variational Monte Carlo calculations with a Gutzwiller wave function and are 
generally somewhat larger (about 5 %) than those obtained in the HFA. The exact 
results are estimated to be smaller than the HFA results. This discrepancy is 
due to the neglect within the SBMF approach (as in the HFA) of the current- 
current correlation part of Ps- The remaining kinetic part of Ps agrees very well 
with Quantum Monte Carlo (QMC) calculations of the kinetic energy or effective 
hopping integral. Also off half-filling, where our expression for ps no longer has the 
interpretation of helicity modulus of the ground state, it is found to represent the 
effective hopping much better than the HFA. 
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A number of extensions of the present work is possible. Most of the extensions 



discussed below have already been performed within the HFA [|rT| and since the 
SBMF approach turns out to be a renormalized form of the HFA, the qualitative 
effect of such expansions on the present SBMF results can be predicted. A possible 
sequel to the present work (which was not attempted before for the HFA) is the 
calculation of the helicity modulus for the spiral or domain-wall phases, which su- 
persede the AF ground state that one has at half-filling and which is the starting 
point for our calculated helicity modulus. Another extension is to introduce a ho- 
mogeneous magnetic field h in the AF phase at half-filling. Using the well-known 
mapping between repulsive and attractive Hubbard models (see e.g. Ref. the 
corresponding expression for ps than equals the superfluid density of the attrac- 
tive {U < 0) Hubbard model off half-filling (without a magnetic field). A final 
straightforward but tedious extension of the present results is to allow for finite 
temperatures. The formalism set up above is perfectly capable of dealing with this 
more general case, but the calculations become somewhat more tedious than for 
T = 0. 
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Appendix A: Band-renormalization factor 

In this appendix expressions for partial derivatives of the band-renormahzation 
factor for the ferromagnet are given. The formula for q^r, which is a function of 
density n, magnetization m, and density of doubly occupied sites (f, is repeated 
here (see (p!^): 



1 - n + <P){n + am- 2<P) + d^/' 



n — am 



2(P 



(n 



am) [1 — i(n + am)] 



(A.l) 

For n, m and d in the physically relevant range (e.g., m should be less or equal than 
n and d"^ should be less or equal than in) q„ attains values between and 1, so the 
free electron bands are narrowed in the SBMF approximation. 

First partial derivatives go-n, Qam and go-d with respect to n, m and d, respectively, 

are: 



qan 
qam 



e 

r„ 



Tg ^ d 
e r_„ 



a 



2d 



e 

r^ r 
r 



d 



(7-^(1 -n - am) 



2e ^ r_o 
e r„ d 



2d 



where we have introduced the abbreviations: 



(A.2) 
(A.3) 
(A.4) 



Ng 
Dg 



1 — 72 + d'^){n + am — 2d'^) + d^/ n — am 
n + am — \{n + amY 



n 



c/2 



\/n + am — 2d'^ 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



To determine where the susceptibility diverges (see Appendix B) , we also require 
qmm, the second derivative of q^ with respect to m calculated at m = (which is 
independent of a): 



1 + 



2(1 - n) 



Dn 



2(1 -n) 

^0 



2ed 



{n - 2d^)Do 



(A.9) 



with 



A^o = Vn - 2^2 (^y/i-n + d"^ + d 
Do = n{l -n/2) . 



(A.IO) 
(A.ll) 
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Appendix B: Homogeneous magnetic susceptibil- 
ity for paramagnet 

In this appendix we derive a formula for the homogeneous magnetic susceptibihty 
X in the paramagnetic phase: 



The derivation proceeds as follows: in the paramagnet, we have m = if /i = 
in the equations (|l7|)-(pl]) of Section 2.1. Then also A = 0. The solutions of the 
consistency equations of the remaining variables we call rio, and Jiq. Now we 
apply an infinitesimal magnetic field 5h. Then m and A will acquire small non-zero 
values 5m and 5\ and the other quantities will deviate slightly from their values for 
= 0, because all are coupled through ([T7|)-(pTD. If we now work at a fixed density 
(i.e., n is not allowed to deviate from hq and equation (^) becomes irrelevant), 
we have four equations containing five small quantities. From these the required 
ratio ^ is obtained. Working out this procedure, by expanding all equations to 
first order in the small quantities, it turns out (perhaps not surprisingly) that the 
equations for 5ji and 5d decouple from those for 5m, 5\ and 5h. Here we give the 
equations for the latter three quantities for the case T = (for which the derivative 
of the Fermi-Dirac distribution is a convenient delta-function): 



5m = a5m + Xo5X , 
5\ = 5h + b5m + a5X , 

where we have introduced the following notation: 

2AfF 



Xo 



with 



2AfFqmP'0 

g3 



ATf = Afino/q) 

fp-o/g 

Sq = I deH{e)e 

J —oo 



B.2) 
B.3) 



B.4) 
B.5) 
B.6) 

B.7) 
B.8) 



and q and qm are the functions q^ and qam taken at m = (see Appendix A). Solving 
and (p.3| ) for x finally gives the required formula for the susceptibility: 



A similar result was given in Ref. ||10||, although there is a factor of 2 difference in 
the SQ-teivd in ( p.6|) . 
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Appendix C: Band-renormalization factor 

In this appendix expressions for partial derivatives of the band-renormahzation 
factor Qs for the antiferromagnet are given. The formula for qg, which is a function 
of density n, staggered magnetization parameter mg, and density of doubly occupied 
sites (P, is repeated here (cf. (p8D-(p^)): 

qg{n, m^, d) = z{n, m^, d)z{n, -m^, d) , (C.l) 

where 



[l-n + d^){n + mg- 2d^) + dy/ 
z{n,ms,d) = 1 = . (C.2) 

+ m,)(l-^ ~ 



n 



Introducing the abbreviations N± and the partial derivative of z with respect to 



iV± = W(n±m,)(l-^^4^), (C.3) 



dz 1 \ 1-n + d^ d z{n,ms,d) i i /r^ ^\ 

— ' ' - [1 — n — rrigj > (C4) 



dm^ 2N+ 1 V n + ms - 2c/2 ~ - 2d'^ N+ 

the first partial derivatives with respect to rris and d are (since we always work at 
fixed density, the derivative with respect to n is not needed): 

a - zin-m ^N ^^(^^^s,rf) , dz{n,-m,,d) 

" N+N.\ Jn2-ni^-4nd^ + Ad^ ^ 2dVl-n + d^ '^^ 



For the calculation of the staggered susceptibility (see Appendix D) the second 
derivative of Qs with respect to mg at mg = is required; the formula is: 

fd^qA An'^-Sn + Sr, -r i2 ^ 2(1 - n + 2^^) 

^= ,3(2-^)3 W^-n + d^ + d\ (-2^-^(2 -n)(n- 2^2)- 

(C.7) 

Appendix D: Staggered magnetic susceptibility for 
the paramagnet 

In this appendix we derive a formula for the staggered magnetic susceptibility in 
the paramagnetic phase Xs- 
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The procedure is as follows: the consistency equations (|3^)-(|36D for the antiferro- 
magnet in Section 2.2 allow for a paramagnetic solution in which hg = leads to 
ms = as well as Ag = 0. Starting from this solution, we apply (at a fixed density) 
an infinitesimal staggered magnetic field Shg. This introduces small changes in the 
other parameters, in particular mg and As acquire small values SrUs and 6Xs (the 
changes in d and fl are irrelevant for the present discussion). Restricting ourselves 
to the T = case, the two equations relating Shg, Srug and 6Xs are (expanding the 
consistency equations to first order in 6hs, 6ms and 6Xs): 

6ms = Xs,o6Xs , (D.2) 
6Xs = 6hs + bs6ms , (D.3) 

where the following abbreviations are introduced: 

2 Af{e) 
Xs,o = - de p.4 
q J-oo e 

bs = 2q^^m, / de^{e)e . (D.5) 



Here the parameter q is qs taken at mg = (see Appendix C). Solving (|D.2|) and 
(p.3|) for Xs one has: 

Sms Xsfl 
6hs 1 - 6sXs,o 
A similar result was obtained in Ref . | T0[ . 
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Tabic I. Comparison of location of tripod points (where three phases meet) in 
{4:t/U, 6) phase diagram between Hartree-Fock approximation (HF) and slave-boson 
mean-field approximation (SBMF). Also the critical hole densities 6c for antifer- 
romagnetism (AF) and ferromagnetism (F) to occur in both approximations are 
compared. The last column gives the ratio of the HF and SBMF result in each 
case. 
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,623639 


0, 


.951682 


8 


2.98950 


0, 


.863161 


0, 


.235612 


0. 


,485104 


0, 


.959478 


10 


4.11613 


0, 


.913012 


0, 


.193926 


0. 


,393528 


0, 


.968953 


12 


5.23030 


0, 


.940510 


0, 


.163654 


0. 


,330002 


0, 


.976425 


16 


7.39950 


0, 


.967337 


0, 


.123965 


0. 


,248782 


0, 


.985641 


20 


9.51220 


0, 


.979394 


0, 


.099524 


0. 


,199420 


0, 


.990505 


200.1 


100 


0, 


.999800 


0, 


.009998 


0. 


,019976 


0, 


.999900 



Table II. Self-consistent parameters for antiferromagnetic ground state at half-filling 
in slave-boson mean-field approximation as a function of U/t. The ground-state 
energy for = equals — 16/7r^ exactly. 
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U/t pH'^A pf pGWVMC 






0.2026 


0.2026 





1 


0.2023 


0.2012 




2 


0.1960 


0.1953 


0.197 


3 


0.1820 


0.1847 


0.185 


4 


0.1650 


0.1713 


0.172 


6 


0.1332 


0.1421 


0.141 


8 


0.1090 


0.1162 


0.117 


10 


0.0912 


0.0962 


0.098 


12 


0.0781 


0.0814 


0.083 


16 


0.0602 


0.0618 




20 


0.0488 


0.0497 





Table III. Spin-stiffness ps of tlie antiferromagnetic ground-state at lialf-filling as 
a function of U/t as calculated in tlie Hartree-Fock approximation (HFA), tfie 
slave-boson mean- field approximation (SBMF), and from variational Monte Carlo 
calculations using the Gutzwiller wavefunction (GWVMC). For U/t = 0, ps equals 
2/7r2 exactly. For HFA and SBMF, results are for an infinitely large lattice, for 
GWVMC, results are for an 8 x 8 lattice, except for U/t = 2,3 which are for 
20 X 20 and 14 x 14 lattices, respectively (see also Ref.|TT[]). 
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Figure captions 

Figure 1 Ground-state {4:t/U, S) phase diagram of the Hubbard model on a square 
lattice, restricted to simple magnetic phases: paramagnet (PM), antiferromagnet 
(AF), weak (WF) and strong (SF) ferromagnet. S is the density of holes: 1 — n. 
(a) Construction diagram showing all continuous and first-order transition lines 
obtained in slave-boson mean-field approximation (SBMF), (b) phase diagram in 
SBMF, and (c) corresponding phase diagram in the Hartree-Fock approximation. 
Note the difference in scales of (b) and (c). 



Figure 2 Helicity modulus ps for the repulsive Hubbard model on a square lattice 
at half-filling as a function of U/t. Shown are results from the Hartree-Fock approx- 
imation (HFA), from the slave-boson mean-field approximation (SBMF), and from 
variational Monte Carlo calculations using a Gutzwiller projected wave function 
(GWVMC, from Ref.|ll). 



Figure 3 Effective hopping integral tes/t for the repulsive Hubbard model on a 
square lattice at half- filling as a function of U/t. Shown are results from the Hartree- 
Fock approximation (HFA), the slave-boson mean- field approximation (SBMF), and 
Quantum Monte Carlo calculations (QMC, from Ref. ||32||). 



Figure 4 Effective hopping integral tes/t for the repulsive Hubbard model on a 
square lattice as a function of electron density n. (a) For U/t = 4, results are shown 
from the Hartree-Fock approximation (HFA), the slave-boson mean- field approx- 
imation (SBMF), and Quantum Monte Carlo calculations (QMC, from Ref.[Q), 
and (b) for U/t = 4, 8, 12, and 16 from SBMF calculations (QMC results for U/t = 4 
only). The dashed lines indicate the electron densities Uc for which in the phase 
diagram in Fig. 1(a) a continuous PM/AF transition takes place. For U/t = 8, 12, 
and 16 the value of Uc is (approximately) equal to 0.8 in each case. 
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